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We study the origin of buoyancy forces acting on a larger particle moving in a 
granular medium subject to horizontal shaking and its corrections before fluidiza- 
tion. In the fluid limit Archimedes' law is verified; before the limit memory effects 
counteract buoyancy, as also found experimentally. The origin of the friction is an 
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CN . excluded volume effect between active particles, which we study more exactly for 
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^ \ a random walker in a random environment. The same excluded volume effect is 
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also responsible for the mutual attraction between bodies moving in the granular 
medium. Our theoretical modeling proceeds via an asymmetric exclusion process, 
i.e., via a dissipative lattice gas dynamics simulating the position degrees of freedom 



<~j ' of a low density granular sea. 
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I. INTRODUCTION 



A main characteristic of granular media is that their behavior varies between being more 
fluid- and being more solid-like. Initializing flow via shaking or stirring fluidizes granular 
baths. That has been observed in various experiments and simulations, in particular by 
verifying Archimedes' law, jl-3]. Also phenomeno logical arguments have been added to the 
understanding of the buoyancy force in granular media, e.g. from using the Enskog hy- 
drodynamic equations Qj. Nevertheless numerous controversies have remained and various 
corrections must be considered. In the present paper we take up a simple excluded vol- 
ume model to study the origin and the corrections to Archimedes' law around the fluid limit. 



Fluid-like behavior of a granular material should obviously include hydrostatic and hy- 
drodynamic effects. A natural way to study the origin of granular hydrostatics and their 
possible corrections is via flow induction, i.e., stirring or shaking the system. That causes 
energy transfer to the grains which is dissipated again in the collisions between the moving 
grains. Further simplifications can help to understand the essence of what happens. In that 
spirit we consider the asymmetric exclusion process to simulate the dynamics of the grains 
(monomers) with one large particle (rod) immersed in them. The condition of detailed 
balance enables us to identify the buoyancy force on the rod, as function of its size and of its 
relative weight, and locally as function of the height. We recover Archimedes' law in the fluid 
limit, and we create a theoretical framework for a detailed study of possible corrections. 
Corrections arise from various effects such as from the discrete nature of the lattice where 
the lattice spacing measures the size of the grains. Corrections also arise from thermal 
effects especially when the grains are themselves immersed in a heat bath (e.g. hot air), and 
as studied here, from finite shaking rates. Other possible corrections arise from convection 
currents in the granular medium, to which we turn briefly in Section IIVI 



To the extent that our modeling via the asymmetric exclusion process is relevant for 
the experimental conditions of granular media under shaking, Archimedes-like behavior was 

n n 

already predicted in [5|. Five years later Archimedes' law was confirmed experimentally [1| 
in a granular medium, and some corrections were explored before the fluid limit. In the 
mean time further experimental work and simulations which show an explicit Archimedes 



like behavior have been added, including [2|, |3[ . We come back to the basic set-up 



A. Phenomenon & Experiments 

The experiment [lj] consists of a bi-dispersion of glass beads which form the granular 
bath. This bath is placed in a rectangular box, on which horizontal shaking is applied 
through vibrations on the walls. Gravity controls the vertical motion. The bath ffuidizes 
for larger energy input through the external shaking. The energy of shaking is quantified 
by a dimensionless quantity T oc f 2 A/g related to frequency / and amplitude A of shak- 
ing which in the case of horizontal shaking exercise a similar effect. One observes a clear 
transition; before some energy threshold is reached the state is mostly solid and there is a 
clear boundary where the unfluid state becomes fluidized and the principle of Archimedes 
is obeyed henceforth. The correction to Archimedes' law can be described by a drag force. 
Several experiments have studied the behavior of this drag force l|, |2|, |6|-|8| . Most agree that 
in the fluidized and low density limit the drag force is observed to be proportional to the 
velocity. Moreover the coefficient of viscosity is seen to vary exponentially with the shaking 
amplitude If]], or, in observed inversely proportional to the shaking amplitude. For 
denser media the drag force was seen to be varying logarithmic with velocity 3]. In our 
model we discuss the nature of drag force and corrections to it in Section IIIIBl 



B. Model and results 



The goal of this paper is to give a simple theoretical model giving a natural interpretation 
of these experimental results, and at the same time to allow further studies on modifications 
of Archimedes' principle under the influence of fluctuations and nonequilibrium effects. On a 
more fundamental level, our model illustrates the construction of statistical forces. We have 
here an example of how the effective motion of the bigger particle gets changed indirectly 
by the influence of that same particle on the smaller particles. In other words, the rod 
(bigger particle) is active in changing its environment. The result can be characterized as 
an additional effect modifying Archimedes' buoyancy before the fluid limit. 

The granular medium that is considered are particle conserving lattice gases that evolve 
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under excluded volume conditions. The configuration is completely determined by the par- 
ticle occupations (there are no velocities) but the dynamics is dissipative and therefore 
simulates to some extent the behavior of low density granular materials. 
The plan of the paper is as follows. In the next section we describe our model, introduced 
in j^l which is an exclusion process of walkers on a regular lattice in the presence of a spa- 
tially extended object. In Section III Bl the fluid limit yields Archimedes' law in terms of a 
conservative buoyancy force proportional to the local fluid density times the volume of the 
object. Further the continuum limit to the spatial degree of freedom reveals the first level of 
corrections to the nature of buoyancy force. Fluctuations due to the surrounding heat bath 
bring in the next level of correction to the usual hydrostatic formulation. 
The main result appears in Section III II where we investigate corrections to the fluid limit. 
We map our model to a simpler one-dimensional problem, to obtain the motion of a ran- 
dom walker in a stochastic dynamical environment. There, the correction to Archimedes' 
law is explained by a memory effect which acts as friction on a rising object. Continuum 
corrections bring additional clarification. We discuss how our model adds a simple heuristic 
picture to the experimental observations of lj. The last sections discuss additional aspects 
that become visible through our theoretical modeling. A discussion on the existing mecha- 
nisms of intruder dynamics in granular media and a comparison with our model appears in 
Section IIVI In section IIVBI results of simulations modeling our system are shown. Future 
possibilities of nonequilibrium effects and collective behavior are alluded upon in the last 
sections. 



II. A ROD IN A LATTICE FLUID 



We consider a stochastic dynamics for the motion of a rod in a lattice fluid composed 
of monomers. The model was introduced in jlj], obtaining a Markovian reduced dynamics 
for the rod motion in the limit where between any two moves of the rod, the monomer 
configuration has the time to relax. That can be interpreted as the fluid limit in which the 
rod always "sees" the fluid of monomers in equilibrium. 
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FIG. 1. Two-dim model 
A. Model 

All motion takes place on the square lattice Z 2 where the mesh size, taken unity here, 
gives the size of the grains. These grains (also called, monomers) can occupy sites % = (x, y) 
having "vertical" coordinate y and "horizontal" coordinate x. There can be at most one 
grain per site. There is one big particle or big grain, called rod and we only follow its vertical 
position. The vertical position of the rod at time t is denoted by Y t taking values in Z. The 
rod occupies N £ {2, 3, . . .} lattice sites in the horizontal direction, i.e., the region 

A N (y) = {(0,y),{l,y),...,{N-l,y)} (II.l) 

is forbidden for the monomers when Y t — y. 

The horizontal jumps of the monomers are symmetric at rate 7. Increasing the rate 7 speeds 
up the monomer dynamics in the horizontal direction (orthogonal to the motion of the rod). 
The vertical jumps of the monomers and the rod are asymmetric; modeling a gravitational 
field. Note that the vertical motion is not speeded up, but of course it is influenced by 7 as 
well. 

More formally, the microscopic dynamics looks as follows; see Fig. 1: The monomer 
configuration is denoted by 77 £ {0, l} z2 ; T)t(i) = means there is no monomer at site i at 
time t and r]t(i) = 1 if there is a monomer at site % at time t. The dynamics is of exclusion- 
type because all motion is via jumping to vacant sites, and rod and monomers never overlap. 
A monomer moves horizontally to a vacant nearest neighbor site, symmetrically with rate 
7. It moves vertically up with rate p and down with rate q. The rod only moves vertically, 
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up with rate a and down with rate b. We choose p/q, a/b < 1 to represent the gravitational 
field, e.g. via p/q = exp(—mg/kT),a/b = exp(—Mg/kT), where M,m denote the mass of 
the rod, respectively of the monomers; kT is a typical unit of thermal energy at temperature 
T which plays little role in what follows, except for allowing fluctuations. The temperature 
could also refer to an additional heat bath that makes contact with the grains. Yet, granular 
media are typically a-thermal in which case we think of Mg, mg » kT. The lattice unit is 
not indicated; it is taken to be one and should be thought of as the size of the grains. All 
that motion gets summarized in the formal generator 

Lf( V , y) = al[ v (i) = 0, V* e A N (y + 1)] [f( V , y+1)- f(v, y)} 

+ b Ifo(i) = 0, Vi e A N (y - 1)] [f(rj, y - 1) - f(rj, y)] + 

E {pv® ( x - ^ + 1)) Wi, ^ + 1) i Mv)\ x [f(v i>{il ' i2+1 \ y) - f(v, y)\ 

i=(u,«2) 

+ q r)(i) (1 - rjik, i 2 - 1)) /[(•!, i 2 - 1) i A N (y)} x [/(t^ 1 ^, y) - f(r), y)]} 

+ 7 E im^A N ( y ) = ds][fW J ,y)-f(v,y)} (ii.2) 

where /[•] is the indicator function of the event in the brackets, giving one or zero depend- 
ing on the event being realized, and rf^ is the grain configuration j] after switching the 
occupations in sites % and j. The last term represents the horizontal shaking in which the 
occupations of (horizontal) nearest neighbor pairs (ij) get exchanged. Observe that there 
is always both horizontal and vertical motion, subject to the exclusion rule, which, besides 
from the shaking, can arise from an extra heat bath in which the grains are moving and 
with which energy can be exchanged. 

We need (III. 2D for writing down the kinetic equations that are all of the form 

j t (f(riuY t )) = (Lf(r ]u Y t )} 

Here and from now on, brackets (•) are with respect to the stochastic dynamics and over the 
following initial conditions. At starting time t = we put Y t= o = so that the rod starts 
from the center of the lattice, but that is really arbitrary. For the initial distribution on 
monomer occupations we take density 

d(x, y) = d{y) = - K fr/ g ) y (II.3) 
v yi yyi l + K(p/q)y y ' 

for parameter k > 0. This formula satisfies pd(y) (1 — d(y+ 1)) = q d{y + 1) (1 — d(y)) which 
is a detailed balance relation for the motion of the grains, see also (HV.ip below. The density 
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varies between zero (at the top) to one (at the bottom). The height where d(y) = 1/2 scales 
like y ~ log ft. The derivative of the density at that height (where the transition is made 
between higher and lower density) is proportional to mg/kT. Therefore, choosing k ~ 1 
and large mg/kT corresponds to a more constant density as in a liquid or as in a granular 
medium in a container filled from the bottom to around y = 0; on the other hand, looking 
at positive y, for k = 1 and for smaller mg/kT corresponds to a gas condition where (1II.3P 
simulates a barometric formula. Low density granular media under heavy shaking would 
also fall in that category which in our modeling is most typical. The density is constant in 
the horizontal direction (x), but always conditioned on having rj(i) = for all i covered by 
the rod, i.e., for all % G A^{Y). 

More precisely, we let Vd denote the product measure on {0, l} z2 with density 

Probes, y) = 1] = d{y) (II.4) 

defined by (111.3p . The conditional probability 

iy° d = M-\r](i) = oyzeA N (o)) (II.5) 

is then the initial distribution on the monomers. The dynamics such as defined above gives 
rise to the Markov process (r] t , Y t ). 



B. Fluid limit 

In the limit 7 t +00 the motion of the rod decouples from the monomer dynamics. Then, 
the reduced dynamics of the rod becomes that of a random walker with rates directly given 
in terms of the equilibrium fluid density: 

the rod moves up y — > y + 1 with new rate a[l — d(y + 1)]^ and goes down y — > y — 1 with 
new rate b[l — d(y — 1)} N . The factors [1 — d(y ± 1)]^ of course express the plausibility 
of having space for the rod to move from height y to y ± 1; there must be a hole of size 
N. Hence, in the limit of excessive horizontal shaking, the rod is doing a continuous time 
random walk on Z with backward generator 

L™f(y) = a[l - d(y + l)] N [f(y + 1) - f(y)) + b[l - d(y - l)] N [f(y - 1) - f(y)), (11.6) 



8 



The density profile d is obtained from (1II.3P - 



We call this limit 7 | oo the fluid limit. The reason of the decoupling is that the 
monomers relax to their stationary reversible density in between any two moves of the rod. 
The resulting motion flll.6j) is itself satisfying the condition of detailed balance for a potential 
V, which can be interpreted as giving rise to a conservative force F given by the logarithmic 
ratio of up versus down rates 

b[l-d(y-l)] N 



F(y) = -V(y) + V(y-l) = -kT In ■ 



-Mg - NkT ln[l 



a[l-%F 

d(y) - d(y- 1] 



d(y) J 

(II.7) 



To go to a continuum description (at least in the vertical direction) we introduce a lattice 
mesh of size e > 0, under which In a/b = —Mge/kT, In p/q = —mge/kT. This imagines that 
the grains are of vertical size e. We compute the force F £ as e 4 0: 

-V(y) + V { y-s) = N1+ <!M, (II . 8) 

e el — d(y) 

where the density d(y) is now on M and, similar to (III. 3}) . verifies 

d'(y) = -d(y)(l-d(y))^ (11.9) 
Hence, (1II.8P becomes F £ — > F for 

F(y) = -Mg + d(y)mgN (11.10) 

which is Archimedes' law for the total upward force on a body of volume iV and mass M 
replacing a weight equal to mgd(y)N of fluid. For short, we call this the buoyancy force. 
The rod will thus move to a height where the fluid density is proportional to 1/N. That 
is the equilibrium position, consistent with Archimedes' characterization of the hydrostatic 
equilibrium position, {9]. The force ( 1II.7P is the correction to the Archimedes' force (IH.lOp . 
due to the finite size of the grains. 



The motion can be studied in the diffusive limit where we also rescale time ~ s 2 . That 
means to take for example 

p -Mge/{2kT) Mge/(2kT) 
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and to expand the generator 

Lf w f(y) = a [l-d(y + e)f [f(y + e) - f(y)} +b[l-d(y- e)f [f(y - e) - f(y)} 
(see (1II.6j) ) in orders of e. The result is that L^ w f(y) — > Lf(y) with 

Lf(y) = (D f)'{y) + x (y) F(y) f'(y) (11.12) 

with 

F{y) — N d(y)mg — Mg, D(y) = [1 - d(y)} N , X {y) = ^ 
We again have made use of (III. 91) in the continuum limit. The result ( 111.12)) is the generator 



of an overdamped diffusion equation with diffusion coefficient D(y). The corresponding 
Langevin equation, in the Ito-sense, is given by 



Vt = xiVt) F(y t ) + D'(y t ) + y/2D(y t ) & (11.13) 

for white noise £f. The force F is exactly the one found in (IH.lOj) . as in Archimedes' law. 
The diffusion D(y) is related to the mobility xiu) y i a the Einstein equation kT x(y) — D(y). 
The term with the derivative D'(y) is due to the Ito-convention. 

The above analysis concludes that the two-dimensional lattice model on a lattice with a 
simple exclusion dynamics, provides a reasonable description of the hydrostatic behavior of 
granular matter in the limiting (fluidized) case. We have seen above how the discreteness of 
the vertical lattice-direction makes a first correction, easily studied for small lattice mesh. 
A second type of correction is due to fluctuations. The rod dynamics generated by (1H.6|) is 



stochastic. The fluctuations of the rod about its equilibrium position can be studied in the 
large N limit with standard deviation around the mean going as 1/ y/~N. For the third major 
type of correction, we study the approach to the fluid limit. The next section is devoted to 
these questions. 

III. RANDOM WALK IN A DYNAMICAL ENVIRONMENT 



To this end we propose a contracted description of the model, coarse-graining it to an 
effective one-dimensional model; see Fig. 2. The idea is as follows. The essential aspect of 
the monomer dynamics as far as regards the rod, is whether there is a hole above or below 
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FIG. 2. Contracted description 

the rod in which it can jump. We therefore summarize all of the monomer configuration 
T)(x, y) by variables a(y),y G Z, that specify whether or not there is a monomer in the region 
A N (y), see (JO): 

<j{y) = 1 if there is no hole at y 

= if there is a hole at y (III. 2) 

More precisely, there is a hole at y if 77(f) = for all % G A^{y). The position of the 
rod is still denoted by Y. We assume as major simplification that the (contracted) system 
(at(y),Y t ),y G Z, t > 0, undergoes a joint Markov process which mimics the original model 
in the following sense. 

The rod moves up y — > y + 1 with rate a if there is a hole at y + 1, i.e., if a(y + 1) = 0. The 
rod moves down y — y y — 1 with rate b if there is a hole at y — 1, i.e., if cr(?/ — 1) = 0. The rod 
never moves to a position y where there is no hole, a(y) = 1. For the monomer dynamics, 
we assume that the &t{y) flip <=^ 1 with different rates depending on y, and depending on 
the position of the rod. More precisely, cr t (y) has rate q(y) for the change 1—^0 and has 
rates p(y) for — > 1 except when Y t = y because then it must remain zero; see Fig. 3. 
So formally the backward generator of our new Markov process is 

Cf(a, Y) = a[l — a(y + 1)] [f(a, y + 1) - f(a, y)} 
+ b[l-a(y-l)}[f(a,y-l)-f(a,y)} 

+ [t 1 - '(*)] p(y)(i - W + *(y) ?(y)] n - /(^ y)] (iii.s) 



11 




FIG. 3. Contracted dynamics 

where a y is the hole-configuration obtained after flipping the occupation at y: u v {y') = cr(y') 
if y 7^ y' and cr y (y') = 1 — a(y) if y = y'. At no point in time could the region occupied 
by the rod be simultaneously occupied by a monomer and vice versa, so o- t (Y t ) = always. 
This is the influence of the active particle (the rod) on the fast degrees of freedom. 



A. Interpretation 

Obviously, the rates p(y), q(y) must be interpreted in terms of the monomer density d(y) 
at y with their dependence on the size N of the rod and on the amount of horizontal shaking 
7. Comparing fjIII.3|) with (III. 21) suggests further interpretations. 



It remains that a/b = exp(— Mg/kBT) where M is the mass of the rod. We can think of 
the monomers as blinking lights, red (cr(y) = 1) for no passage of the rod, and green (cr(y) = 
0) for passage allowed. In the original two-dimensional model a hole (green light) a(y) = at 
y represents the fact that there are no monomers at the sites A^{y) = {(0, y), (1, y), . . . , (N— 
1,1/)}, and (red light) a(y) = 1 means that some monomer can be found in the region A^{y). 
Abbreviating 

p{y) = 1 - (1 - d(y)f 

we take therefore 

p(y) = ip(y), q(y) = t(! - p(v)) ( IIL4 ) 

for the rate at which a hole gets removed, respectively created. Each depend on N and on y 
but observe that p(y) + q(y) = 7 which is the horizontal shaking rate. The stationary hole 
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density at y for that two-state Markov process becomes 1 — p(y) = (1 — d(y)) N which is the 
correct hole probability in the original monomer-model, cf. (III.4p . Of course the weight of 
the monomer is represented in the density d(y) via fill. 31) . 



Here also, for our simplified model we can take the fluid limit 7 j" 00. By simpler 
arguments than in jfjj], the motion of Y t decouples from that of the at and by our choice 
( 1III.4|) we find exactly the same limiting motion of the rod as given in (111.61) . The two 
models mathematically agree in the fluid model but our second model allows more easily to 
find the most significant contribution before the fluid limit, to which we turn next. 



B. Before the fluid limit 



Suppose we find the rod at time u in position y. Thus, at that time the hole probability at 
y is equal to 1. When the rod jumps from position y to say position y + 1 at time u it leaves 
a "hole" at position y which remains a hole until it gets occupied by either monomers or 
by the rod again. In the fluid limit the monomer dynamics is fast enough and they relax to 
their equilibrium configuration so that at its next jump the rod sees a hole with probability 
1 — p(y). On the other hand if the time of monomer relaxation is longer in comparison to 
the rod, at the next jump at time u + t the rod sees "no hole" with probability 

/0 (y)(l-exp(- 7 t)) (III.5) 

as follows from a simple calculation for the two-state Markov process at(y). Note however 
that the transient density (IIII.5j) is lower than < p(y) and that flIII.5j) is only valid under 



the condition that the rod has not (re-)entered position y during [u, u + t]. It implies that 
at time u + t the rod still sees the hole it left behind at y. We conclude that before the fluid 
limit, the jump rates of the rod depend also on the past and the rod dynamics is by itself 
non-Mar kovian for finite 7. 



We make the above statements now more concrete, and more precise. Knowing the rod's 
motion follows from the evolution equation, for all functions g on Z, 

j t (g(Y t )} = (Cg(at,Y t )} 
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with, from fllP} . 

£<7 (a, y) = a[l - a(y + l)][g(y + 1) - (/(j/)] + 6[1 - a(y - l)][g(y - 1) - g(y)] (III.6) 
Consider for example the expectation 

([1 - a t (Y t + l)]g{Y t + 1)) = ((1 - a t (Y t + 1)\Y„ 0<s<t) g(Y t + 1)) 

where 

(a t (Y t + l)\Y s ,0<s<t) (III.7) 

is the conditional probability of having no hole just above the rod, given the full history of 
the walker (Y s ,0 < s <t). Obviously, history matters. 

Suppose for example that at time t we have Y t = y and that the rod has been there already 
for a time ti. The previous position was either y + 1 or y — 1, from which the rod has moved 
at time t — t\. Before that, at jump time t—ti — ti the rod has been jumping either from y — 2, 
from y or from y + 2, et cetera. In this way the whole history of the rod can be parameterized 
in terms of waiting times and successive positions. We denote such a rod-history by u. Yet, 
the only thing that matters for the expected hole probability at y + 1 in (1III.7P is the last 
time t(oj, y + 1) it was occupied by the rod, since 

([1 - a t (Y t + 1)] | cj) = (1 - d(y + 1))" + p(y + 1) exp(- 7 (t - t(cu, y + 1))) (III.8) 

We put t(uj, y) = — oo if the rod has never been in position y, to realize the initial condition 

(ULU). 

We must now estimate the conditional expectation 

( e -7[*-*(",2/+i)] \(Y s ,se [0, t])=u) (III.9) 

for a history u in which Y t = y. Clearly, that equals exp —yti if before y the rod was at 
y + 1; otherwise (if before the rod was at y — 1) (1III.9I) is certainly less than exp —7^1 + £2), 
which is much smaller than e -7 * 1 for large 7. There are then two cases depending on the 
sign of the rod's "velocity" 

V t = Y t -Y t _ tl (III. 10) 
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We therefore approach the fluid limit by putting ( IIII.9P equal to zero if V t is positive, and 
by putting it equal to 

u(j) = / e-^ 1 (a + b) e - {a+b)ti dt x = — ± (III.ll) 

Jo a + 6 + 7 



if Vt is negative. The integral (IHI.llj) takes the expectation over the exponential waiting 



time distribution for t\. In ( IIII.lip . the sum a + b = v is a good estimate for the average 
speed of the rod, or ^(7) V t = = ^(7) v t , where 1/(7) = l/(v+y) is the friction coefficient 
and v t is the velocity of the rod before it arrived at Y t . 



The drag force on a particle immersed in granular matter was studied in various experiments 
— m a linear dependence on the particle velocity such as proven above was 

observed and corresponds to low density. Inserting this ^(7) we have obtained for large but 
finite 7 that 

j t {g(Y t )) = {L™g(Yt)) 

+ a (p(Y t + l)M7)^%m + 1) " 9(Yt)} 

+ b ( P (Y t - i)M7)^%m - 1) - 9{Y t )\ ) (111.12) 



always with p(y) = 1 — [1 — d(y)\ an d the first line of (1III.12P corresponds to the fluid limit 
(HH. The V t is ±1 as defined in fllTTO]) . 



To recapitulate, the approximation in which we replace f lIII.9p by ( IIII.lip is the following. 



The rate of the rod's dynamics in comparison to the monomer dynamics is such that for the 
rod making a jump y — > y + 1 it can still "see" the gap it left at y + 1 when indeed the rod 
was &ty + l before it came to y. However the rod does not "see" any gaps which were left at 
y + 1 from earlier visits there: the monomer dynamics is fast so that it can erase the trace 
of the rod's trajectory up to one time step ago. The introduction of the velocity V t of the 
rod is a way to re-install the Markov property, where the state of the rod is now defined as 
its position plus its (previous) velocity. In other words, due to the active nature of the rod 
it acquires memory before the fluid limit (which results in a drag force, see below), which is 
most efficiently dealt with by introducing a velocity. 
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For the position of the rod, g(y) = y, 

j t (Y t ) = (a(l - p(Y t + 1)) - 6(1 - P (Y t - 1))) 

+ a (p{Y t + 1)M7)^^) - b (p(Y t - 1)M7)^-^) 

This equation gives the speed of the rod at time t given its current position and previous 
direction V t . If the rod was moving upwards (V t = +1), then it continues moving up with a 
rate a ((1 — p{Y t + 1))) and goes down with a rate (6(1 — p(Y t — 1))) + 6/2(7) (p(Y t — 1)). On 
the other hand if the rod was moving downwards, Vt = —1, then it continues moving down 
with a rate b ((1 — p(Y t — 1))) and goes up with a rate (a(l — p(Y t + l)))+a p{^) (p(Y t + 1)}. 
In comparison to the fluid-limit there is an increase in the rate of return. The rod has a 
higher tendency to go back to the site it started from when the bath is not completely 
fluid. The rate to go forward remains the same as in the fluid limit. This phenomenon of 
a greater tendency to return with possible subsequent oscillations can be interpreted as a 
greater dynamical activity which becomes effective as the bath becomes less and less fluid. 
When the rod has an overall tendency of rising because of the greater buoyancy, the result 
is friction acting downwards. 



As for the fluid limit here also we can make a small mesh analysis and take the diffusive 
limit. We also need to rescale the shaking 7 — > 7/e 2 so that with the choice of (IIII.4I) . 
M7) — 2/(2 + 7) and V t = ±e. We only need to worry about the additional last two lines 
in (1III.12|) . i.e., corresponding to 



e -, (l _M a e 



2kT' 
Mge , 
2kT' 



l-(l-d(y + e)) 
\-{l-d{y-e)) 



N 



N 



2 + 7 
e + 



[g{y + e) - g(y)} 
[g(y - e) - g(y)] 



(ni.i3) 



2 + 7' 

where e~ = 2e if the rod was going down and e~ = when the rod was going up; similarly 
e + = 2e if the rod was going up and e + = when the rod was going down. Again making 
the e— expansion we find, similar to f ill. 12j) . the corrected Langevin equation 

2(1-^)) 



yt = x(yt) F(y t 



with memory term in the friction, v t 



7 



-v t 



D'(y t ) + V / 2£(^) 6 



(111.14) 



yt-yt-dt 



j being the direction of the velocity just 

\\yt~y t -dt\\ & J J 

before time t; the rest of the Langevin equation (1111.140 is interpreted in the Ito-sense with, 
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in particular the left-hand side referring to yt+dt — Ut- The diffusion coefficient remains the 
same as before in the fluid limit, see All. 13[) . The friction ~ (1 — D(y)) increases with higher 
density. 



If we look at the origins of drag or friction in common phenomena like Brownian motion, 
it arises due to a resistance to motion in the form of collisions from the front. The faster 
a tracer particle moves in a thermal bath the more traffic it finds ahead of itself than 
behind. Of course friction appears in all directions against motion and exists at shaking 
of all strengths. The drag force we see here is a variation of this effect. Our system is 
overdamped and nothing of impact or momentum transfer can be discussed; yet interaction 
via excluded volume will be sufficient to generate (another) force which opposes motion of 
the rod. This force appears when granular baths are not completely fluidized. 

The higher the intruder dynamics rate, the stronger its memory of its previous position 
and the greater is the chance of jumping to its old position. The force become weaker and 
weaker in the fluid limit since the memory of the rod is "instantaneously" being wiped away 
by the monomers. We believe that the "drag force" dependence on 7 and intruder velocity 
as seen m away from the fluid limit are explained by this new kind of opposing force 
rather than the conventional understanding of friction in fluids, especially in a low density 
environment where momentum transfer does not play such a big role. 



IV. FURTHER REMARKS 

A. Segregation effect 

Effects of buoyancy in granular media have been widely studied both theoretically and 
experimentally. It is not always easy to distinguish between anti-gravity effects and buoyancy 
as in Archimedes' law. That connects with the variety of segregation effects in granular media 
upon shaking them. Reference 



10[ discusses various mechanisms that can work together for 
the segregation of grains. 

Buoyancy, i.e., rising/sinking due to pressure gradients, dominates when the fluidization 
occurs with no convection. Our model does not show boundary effects and inertia is absent. 
The limiting motion is overdamped in the diffusion limit. 
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Buoyancy is indeed most visible in a vibro-fluidized regime, where only binary collisions 
are prevalent and there is no long time contact between particles. The medium must have 
minimal convection, so the boundaries must be far and the interactions with the boundary 
reduced. In a fluidized regime the effects due to convection as well as inertia are reduced 



enough for buoyancy to be visible, 



10 



llj . In the unfluidized regime, effects like inertia, 



void filling models (true for vertical shaking) and convection are more important. Another 
difference to be noted is that buoyancy is not just a phenomenon of the larger particle 
climbing up to the top of the pile but refers to a specific dependence of height on the 
relative sizes and densities. 



On the other 
nut effect 



rand, the rising of larger objects in a sea of smaller grains due to the Brazil 



11 



13| . arises in several forms and many competing mechanisms influence the 



motion of the larger particles within a bath (shape, size, forces between particles, shaking 
amplitude and direction, interstitial air and humidity). In Q, [jjl a similar model to ours 
was used to investigate the Brazil nut effect. Sometimes segregation is a result of entropic 
forces which are strong in a gravity free regime and when the frictional forces among bath 
particles and between bath and intruder are such that the entropy of a segregated state is 



higher than a highly mixed state, 



16]. 




FIG. 4. Buoyancy force vs Rod length N at 7 = 15, at height y = 180 and mass of bath particles 
m = 0.02, kT and g taken to be unity 
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B. Simulation results 



Since the one-dimensional reduction of the full lattice fluid model has undergone some 
further simplifications, we have tested numerically whether our approximations appear rea- 
sonable. In other words, we have compared the trajectories of the rod in our approximations 
with those of the true model. The simulation was run on a chain of 300 sites and gravity 
g and kT are taken as unity. In this way we could also numerically verify Archimedes' law 
in the large 7 limit. In the fluid limit the buoyancy force varies linearly with the size iV of 
the rod, see Fig. 4. The straight line indicates that for such a large 7 the monomer bath is 
fluid-like. For a given height y buoyancy force is calculated by estimating the weight of the 
rod which would exactly balance the force from the bath at that height. In the fluid limit 
at equilibrium then, the weight of the rod is equal to the upward force(called buoyancy). 

Before the fluid limit is reached buoyancy force varies with 7; for large 7 the buoyancy 
force tends to a steady value as given by Archimedes' law, see Fig. 5. That must be compared 
with Fig.3 in The buoyancy force grows with 7 and after a certain critical value which in 
the simulation was 7 = 4.0, it saturates. Fig. 5 shows the Archimedes' force (III. 101) corrected 



with the friction term as it acts in flIII.141) . with the 2/(2 + 7) kind of variation. 



C. Longer memory 

Instead of considering memory only until one time step before, one could also take two, 
three or more time steps long memory. That means, to consider again (IIII.9I) and to take into 
account contributions from alternative histories. These contributions are all of smaller order, 
with each correction falling as an inverse power of 7. The power arises from performing the 
integral like in (IHI.lip but now the time is a sum of exponential variables, so that we get 
corrections like f^( , y) n - 



D. Nonequilibrium seas 

Another type of correction to Archimedes' law comes from the possible nonequilibrium 
nature of the medium, even in the fluid limit. Indeed, an essential aspect of our analysis 
above was that the sea of grains (monomers) reach their equilibrium between two moves 
of the rod (in the fluid limit), and that equilibrium is the same with or without the rod. 
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0.032 - 




FIG. 5. Buoyancy force vs 7 for N = 20 at height y = 180 and mass of bath particles m = 0.02, 
kT and g taken to be unity 



In other words, the stationary distribution of the constrained dynamics of the monomers 
given the rod's position gives exactly the same as conditioning the stationary distribution 
of the joint monomer-rod dynamics on the position of the rod. That is only valid under 
the condition of detailed balance, see e.g. Lemma 3.1 in In the present paper detailed 
balance is forced by the specific choice of density profile (111.31) for which holds that 

pd(y) (1 - d(y + 1)) = qd(y + 1) (1 - d(y)) (IV. 1) 

In the small mesh (continuum) limit this detailed balance condition ( II V. 1 j) becomes d'(y) = 
—mgd(y) (1 — d(y))/(kT) as repeatedly used in the derivation of the Langevin equations 
( III.13j) and ( IIII.14j) . and in the validity of the corresponding Einstein relation. When detailed 
balance is violated and the granular sea shows an irreversible steady behavior, the motion 
remains much less understood. 



E. Collective effects 



A final source of corrections to Archimedes' law is due to the interaction with other 
intruders (rods). Here we are really speaking about a whole new range of phenomena in 

can occur. 



which pairing of particles, and more general collective effects as flocking 
The main underlying reason however is already visible from the simple analysis of the present 
paper. The excluded volume effect of one rod not only creates a bias for itself to return to 
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the place it was before (creating drag, [19|), but also creates space where another intruder 
can hop into, and thus "attracts" other rods and intruders. This granular-hydrodynamic 
interaction is long range and is expected to be proportional to the speeds of the rods, 
but more exploration is needed. This interaction qualitatively resembles the long range 
hydrodynamic interaction between colloidal particles in suspensions. These interactions 
come through due to the Stokes-like force applied by the suspension on a moving colloid 
which results in long-range interaction between two colloids connected through the Oseen 
tensor. Here, the collective behavior of multiple intruders results from simple exclusion and 
memory tracks left on the bath. 



V. CONCLUSIONS 

There are few "mesoscopic" toy-models of granular effects, and there is a wealth of ex- 
perimental data. We have presented a theoretical model comprised of a large bath of small 
particles and the motion of a large body through this bath. We only consider positional 
degrees of freedom with a dissipative dynamics. Fluidization is controlled by the degree 7 of 
horizontal shaking; it modifies the relaxation time of the bath degrees of freedom. The fast 
and the slow degrees of freedom get separated and the extended body undergoes buoyancy 
according to Archimedes' law. Corrections before the fluid limit can be modeled in terms of 
memory effects, where the big particle is biased to fall back in the hole it left behind. This 
creates friction proportional to the velocity with a coefficient that is inversely proportional 
to the shaking amplitude 

Granular medium is fertile ground for study of origin and behavior of statistical forces 
which arise due to coarse- graining. Two of those are studied here, buoyancy due to the 
pressure gradients and additional memory effects creating friction. As an outlook, the study 
of effective forces and interactions out-of-equilibrium remains very much open. 
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